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ABSTRACT 

We examine in detail the recent proposal that extreme Cosmic-Ray-Dominated-Regions (CR- 
DRs) characterize the ISM of galaxies during events of high-density star formation, funda- 
mentally altering its initial conditions (Papadopoulos 2010). Solving the coupled chemical 
and thermal state equations for dense UV-shielded gas reveals that the large cosmic ray en- 
ergy densities in such systems (Ucr~ fewx(10 3 -10 4 ) UcR,Gai) will indeed raise the mini- 
mum temperature of this phase (where the initial conditions of star formation are set) from 
~10K (as in the Milky Way) to ^50-100 K. Moreover in such extreme CRDRs the gas tem- 
perature remains fully decoupled from that of the dust, with Tkin^Tdust, even at high densi- 
ties (n(H 2 )~4 5 -10 6 cm~ 3 ), quite unlike CRDRs in the Milky Way where T k — T dust when 
n(H2)>10 5 cm~ 3 . These dramatically different star formation initial conditions will: a) boost 
the Jeans mass of UV-shielded gas regions by factors of ^10-100 with respect to those in 
quiescent or less extreme star forming systems, and b) "erase" the so-called inflection point 
of the effective equation of state (EOS) of molecular gas. Both these effects occur across the 
entire density range of typical molecular clouds, and may represent a new paradigm for all 
high-density star formation in the Universe, with cosmic rays as the key driving mechanism, 
operating efficiently even in the high dust extinction environments of compact extreme star- 
bursts. The characteristic mass of young stars will be boosted as a result, naturally yielding a 
top-heavy stellar initial mass function (IMF) and a bimodal star formation mode (with the oc- 
curence of extreme CRDRs setting the branching point). Such CRDRs will be present in local 
Ultra Luminous Infrared Galaxies (ULIRGs) and merger-driven gas-rich starbursts across the 
Universe where large amounts of molecular gas rapidly dissipate towards compact disk con- 
figurations where they fuel intense starbursts. In hierarchical galaxy formation models CR- 
controlled SF initial conditions lend a physical basis for the currently postulated bimodal IMF 
in merger/starburst versus quiescent/disk star forming environments, while naturally making 
the integrated galactic IMFs (IGIMFs) a function of the SF history of galaxies. 

Key words: cosmic rays - dust, extinction - galaxies: starburst - galaxies: star formation - 
galaxies: interaction - ISM: molecules - ISM: supernovae remants. 



1 INTRODUCTION 

Much of the stellar mass in the Universe forms in starbursts (e.g. 
Blain et al. 1999; Genzel et al. 2001; Smail et al. 2002), spectac- 
ular events during which the star formation rate (SFR) of a galaxy 
rises from few solar masses per year, typical of spirals such as the 
Milky-Way (Mckee & Williams 1997), to several hundred solar 
masses per year as in the local Ultra Luminous Infrared Galax- 
ies (ULIRGs) (e.g. Sanders & Mirabel 1996; Genzel et al. 1998; 
Sanders & Ishida 2004) where merger-driven starbursts take place 
in very compact (D~ 100-300 pc) dense gas disks (e.g. Downes & 
Solomon 1998; Sakamoto et al. 2008). In the so called sub-mm 



galaxies (SMGs), similarly merger-driven starbursts in the distant 
Universe, SFRs can even reach thousands of solar masses per year 
(Smail et al. 1997; Hughes et al. 1998; Eales et al. 1999). While 
stars form invariably out of molecular gas throughout the Universe 
(e.g. Solomon et al. 1997; Frayer et al. 1998, 1999; Greve et al. 
2005, Tacconi et al. 2006), it is in its densest, UV-shielded, and 
Cosmic-Ray-Dominated phase where the star formation initial con- 
ditions, and the characteristic mass of the emergent young stars 
M^ are truly set (Bergin & Tafalla 2007; Elmegreen et al. 2008). 
Indeed, while Photon-Dominated-Regions (PDRs) contain the bulk 
of the molecular gas mass in galaxies (Hollenbach & Tielens 1999; 
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Ossenkopf et al. 2007), their widely varying physical conditions 
remain always far removed from those prevailing in UV-shielded 
CR-dominated dense gas cores. The location of the latter deep in- 
side molecular clouds, a strong temperature regulation via the on- 
set of gas-dust thermal equillibrium at densities n(H2)>10 4 cm -3 , 
and the nearly complete dissipation of supersonic turbulence (Lar- 
son 2005; Bergin & Tafalla 2007; Jappsen et al. 2005), all help to 
keep the thermal state of typical CRDRs a near invariant over a 
wide range of ISM environments. 

Cosmic rays (CRs), accelerated in shocks around supernovae 
remnants (SNRs) or winds associated with O, B star clusters (e.g. 
Binns et al. 2008), are responsible for the heating, ionization level, 
and chemical state of such dense UV-shielded regions in molecu- 
lar clouds (e.g. Goldmith & Langer 1978; Goldsmith 2001). More- 
over, while their heating contribution relative to that from the (far- 
UV)-induced photoelectric effect in PDRs, or turbulence, can vary 
within a galaxy and individual molecular clouds, CRs remain re- 
sponsible for setting the minimum temperature attainable in the 
gaseous ISM (e.g. Goldmith & Langer 1978; Bergin & Tafalla 
2007; Elmegreen et al. 2008). For the Galaxy T k (min)~10K, 
reached in UV-shielded cores with n(H 2 )~(10 4 -10 6 ) cm" 3 . 

A recent study has shown that the average CR energy den- 
sities (Ucr) in compact starbursts typical in ULIRGs (recently 
found also in an SMG at z~2.3, Swinbank et al. 2010), will 
be boosted by the tremendeous factors of UcR~(few)x(10 3 - 
10 4 ) UcR.Gai, potentially transforming their ISM into extreme 
Cosmic-Ray-Dominated Regions (CRDRs). This will dramatically 
alter the thermal and ionization state of UV-shielded cores, and 
thus the initial conditions of star-formation, throughout their con- 
siderable molecular gas reservoirs (~10 9 -10 10 Mq). Higher char- 
acteristic mass M'^' for the young stars formed in such ISM en- 
vironments, and longer ambipolar diffusion timescales for the CR- 
heated UV-shieled cores could thus be a new paradigm for high- 
density star formation in the Universe (see Papadopoulos 2010 for 
details). It is worth mentioning that recently such large Ucr boosts 
(~ 10 3 ' 5 x Ucr, Gai) have been found in the Orion A bar, and while 
not representing the average SF initial conditions in the Galaxy, 
they do allow local studies of the effects of large CR energy den- 
sities on the ISM (Pellegrini et al 2009). Finally CRs have been 
shown to effect the thermal state of primordial gas, catalyzing H2 
formation in the absence of dust ( Jasche, Ciardi, & Ensslin 2007). 

In this work we solve the coupled chemical and thermal state 
equations for UV-shielded gas and its concomitant dust in extreme 
CRDRs without the approximations employed by Papadopoulos 
(2010). We then use the thermal states computed for such a gas 
phase to examine the effects on M'j^' and the so-called inflection 
point of the Effective Equation of State (EOS), set at gas densities 
above which Tkm~Tdust (e.g. Spaans & Silk 2000; Larson 2005 
and references therein). Both of these play a crucial role in deter- 
mining the stellar initial mass function (IMF) in galaxies, which 
could thus turn out to be very different in the CRDRs of extreme 
starbursts. Throughout this work we make the standard assump- 
tion of UcROCp s f r (SFR density) as in Papadopoulos 2010, though 
we now discuss possible caveats (section 4.4). Finally, the different 
IMFs expected in extreme CRDRs will modify the total SFR and 
p s f r values computed in galaxies (from e.g. cm non-thermal or IR 
thermal continuum). In this new context and throughout this work 
SFR and p s f r will now denote these quantities only for the mas- 
sive stars (which define Ucr and power the IR/cm continuum of 
the ISM), with their total values necessarily remaining uncertain, 
depending on the exact IMF, which will no longer be an invariant. 
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Figure 1. The temperature of UV-shielded dense gas cores computed from 
Equation 2 for the range of densities typical in molecular clouds. The 
CR energy densities range from the quiescent environment of the Galaxy 
(U C R=(0.2-l)xU C R,Gal) to those in extreme CRDRs (U C r~(10 3 - 
10 4 )xUcr Gal. Papadopoulos 2010). A dust temperature T(j ust =8 K and 
a Galactic CR ionization rate of £cr = 5x 10~ 17 s" 1 have been assumed. 

2 THE MINIMUM GAS TEMPERATURE IN CRDRs 

In this section we compute the equillibrium gas temperatures for 
the UV-shielded regions of molecular clouds immersed in the large 
cosmic ray energy densities of the extreme CRDRs expected in star- 
bursts with high SFR densities. 

2.1 A simple method 

We first consider the simple case where the chemical state of the 
UV-shielded dense gas cores remains unaltered and the only ef- 
fect of large Ucr values is to heat them much more than in quies- 
cent ISM environments. Following standard practice (e.g. Goldmith 
2001) the temperature of such regions can be estimated from the 
thermal balance equation 

r C R = A Une + A gd , (i) 

where Fcr oc (cr n(H2) is the CR heating with CcrocUcr being 
the cosmic ray ionization rate. The cooling term Aii nc is dominated 
by CO and other molecular lines, while A g d accounts for gas cool- 
ing due to the gas-dust interaction. 

A significant point of departure from the previous analysis by 
Papadopoulos (2010) is the inclusion of the atomic fine structure 
lines OI (at 63 pm) and CII (at 158 pm) in Aii ne . Indeed for the 
CR-boosted temperatures expected in the UV-shielded dense gas 
of extreme CRDRs, and unlike the corresponding gas phase in the 
quiescent ISM of the Galaxy, the cooling from the fine structure 
atomic line of OI (63pm) can also become important and along 
with the gas-dust interaction overtake the CO lines for typical core 
densities of n(H 2 )^10 5 cm -3 when T k >50K . Finally ISM-CR 
interactions can generate internal UV radiation which destroys CO, 
producing traces of CII, thus cooling via its fine structure line at 
158 pm must also be included. The thermal balance equation en- 
compassing all physical processes that can be important in the UV- 
shielded ISM of extreme CRDRs thus is: 
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Aco(T k ) + A gd (T k , T dust ) + Aoi6 3 (T k ) + A C n(T k ) = Tor, 

(for the expressions used see the Appendix). The solutions are 
shown in Figure 1 from where it becomes obvious that the mini- 
mum temperatures attainable in dense gas cores in the ISM of such 
starbursts are much higher than in the Galaxy or galaxies with mod- 
erate p s f r values. Turbulent gas heating was not considered since we 
want to compute the minimum temperature values for UV-shielded 
gas and since supersonic turbulence (and thus shock-heating) has 
usually dissipated in dense pre-stellar gas cores where the initial 
conditions of SF are set. 

From Figure 1 it can be seen that the minimum gas tempera- 
tures possible in extreme CRDRs are lower than those computed by 
Papadopoulos (2010), a result of the extra cooling power emanating 
from the now included atomic transitions that become excited at the 
high temperatures now possible in an ISM permeated by large CR 
energy densities. Nevertheless these minimum possible gas temper- 
atures remain ~5-10 times higher in extreme CRDRs than those in 
the Milky Way or galaxies with moderate p s f r values. Any IR ra- 
diation "leaking" deep in the cloud and warming the concomitant 
dust above ~(8-10) K, or any remnant turbulent gas heating (Pan & 
Padoan 2010) can only raise these minimum temperatures further. 

2.2 Solving for the coupled chemical and thermal state of 
UV-shielded gas cores in extreme CRDRs 

The thermal balance equation for UV-shielded, CR-heated gas 
(Equation 2) and its solutions (Figure 1), while physically trans- 
parent, they do not account for the fact that ISM chemistry (and 
thus the abundance of coolants) is also strongly CR-regulated in 
such regions. This makes the problems of thermal and chemical 
balance of CR-dominated gas strongly coupled that must be solved 
together, especially given the very large CR energy densities pos- 
sible in the CRDRs of compact starbursts. We do so by utilizing 
the advanced time-dependant UCL-PDR code (the latest version 
described in Bell et al. 2007) to compute gas temperatures of CR- 
heated molecular cloud interiors self-consistently with their CR- 
regulated chemistry for a wide range of CR energy densities. This 
model has already been successfully used several times to model 
PDRs in the Galaxy (e.g. Thi et al. 2009). 

The heating mechanisms included are: (i) photoelectric ef- 
fect from silicate grains and Polycyclic Aromatic Hydrocarbons 
(PAHs), (ii) H2 formation on gain surfaces, (iii) H2 photodisso- 
ciation, (iv) H2 UV fluorescence, (v) CII recombination, (vi) in- 
teraction of low-energy cosmic rays with the gas. The latter dom- 
inate gas heating deep inside H2 clouds (A v ^5) where dense gas 
cores reside and the initial conditions for star formation are set. At 
cloud surfaces hot gas emits Lya and OI 63/im lines while deeper 
in (A v >0.1), gas cooling occurs via fine-structure line emission 
of oxygen and neutral and singly ionized carbon. Vibrational and 
ro- vibrational lines of H2 and CO also cool the gas while CO rota- 
tional lines are excellent coolants for the cold inner regions as the 
first rotationally-excited level reaches down to 5.5 K. The interac- 
tion between the gas and the dust grains acts as cooling or heating 
agent for the gas depending on the difference between the gas and 
dust temperature. H2 and CO self-shielding are taken into account. 
The code calculates the gas chemical abundances, emergent cool- 
ing line fluxes, and temperatures self-consistently at each depth in 
the cloud, and it has been benchmarked against other state-of-the 
art photodissociation region codes in Rollig, M. et al. (2007). 

We set a nominal incident far-UV field of G D = 1 in Habing 



units (expected to be even lower for the dense pre-stellar cores 
deep inside H2 clouds) and a residual micro-turbulent line width of 
AV m _ t urb=5kms _1 (which enters the expressions of line optical 
depths). In all cases chemical and thermal equillibrium is reached 
fast, within T c hem,th^ 10 4 yrs, and is always shorter than the lo- 
cal dynamical timescale which drives the ambient star formation 
events and thus the temporal evolution of Ucr- The results, shown 
in Figure 2, corroborate the high gas temperatures obtained using 
the simpler method. This reveals the cosmic rays as the most im- 
portant feedback aspect of star formation capable of significantly 
altering its initial conditions during events of high star formation 
density via the strong heating of the deeply embedded and UV- 
shielded dense gas regions of molecular clouds. It can thus oper- 
ate unhindered in the high dust extinction environments of extreme 
starbursts such as local ULIRGs or high redshift SMGs. 



2.2.1 Testing the results for stronger radiation fields 

We also ran the same grid of models with the incident far-UV ra- 
diation field scaling exactly as the CR energy densities which cor- 
responds to the maximum incident radiation field that can be ex- 
pected for dense gas cores (in practice pre-stellar cores are always 
UV-shielded deep inside molecular clouds). The dust temperature is 
now more important than in the previous case of Go=l, and is com- 
puted as part of the code. Its value inside the cloud can influence: 
the cooling by emitting IR photons which interact with line radia- 
tive transfer, the H2 formation by changing the sticking probability, 
the evaporation/condensation of molecules, and most importantly 
the gas-grain heating/cooling. The UCL-PDR code uses a modified 
version of the formula of Hollenbach, Takahashi, & Tielens (1991) 
(hereafter HTT) corrected for the various mean grain sizes. 

The HTT formula has been further modified to include atten- 
uation of the far-IR radiation. The incident FUV radiation is ab- 
sorbed and re-emitted in the infrared by dust at the surface of the 
cloud (up to A v ~l mag). In the HTT derivation, this FIR radia- 
tion then serves as a second heat source for dust deeper into the 
cloud. However in their model this secondary radiation component 
is not attenuated with distance into the cloud, remaining undimin- 
ished with depth and leading to higher dust temperatures deep into 
the cloud, which in turn heats the gas to unrealistic temperatures. 
Attenuation of the FIR radiation has been introduced by using an 
approximation for the infrared-only dust temperature from Rowan- 
Robinson (1980), his Equation 30b: T=To(r/r o ) _0 ' 4 where r D is 
the cloud depth corresponding to A v =l mag outer 'slab' that pro- 
cesses the incident FUV radiation and then re-emits it in the FIR. 

From Figure 3 it can be readily seen that the gas tempera- 
tures remain invariant beyond A v ~5 even with with incident far- 
UV fields as enhanced as the CR energy densities. This essentially 
leaves CRs as the only strong feedback factor of star formation op- 
erating efficiently deep inside molecular clouds. Moreover, intense 
CR heating of UV-shielded gas strongly decouples gas and dust 
temperature even at high densities (Figures 1 , 2, 3) with important 
implications on the gas fragmentation properties and thus on the 
emergent stellar IMF (section 3.2). 



3 NEW INITIAL CONDITIONS FOR STAR FORMATION, 
AND THE STELLAR IMF IN EXTREME CRDRs 

The new detailed thermo-chemical calculations presented here sup- 
port the results from the previous study of CRDRs and the effects 
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Figure 2. The gas temperature profiles versus optical extinction inside molecular clouds immersed in a range of CR energy densities, estimated using the 
time-dependant UCL-PDR code solving both for thermal and chemical equillibrium at each point of the cloud. The CR energy densities range from Galactic 
values ~(0.2-1)xUcr Gal. to low-level star forming systems ~10xUcr Gal. and extreme starbursts >~(10 3 -10 4 )xUcr Gal- The densities encompass the 
range typical for star-forming molecular clouds, and the incident far-UV field assumed is Go=l (in Habing units). The dust temperature T(j ust (A v ) profile 
is indicated by the dashed black line. Stronger radiation fields reaching the dense pre-stellar cores deep inside molecular clouds (e.g. because of "leaking" 
far-UV radiation from 0,B stars inside the typically clumpy molecular ISM), can further raise the gas temperatures for A v ^5, while deeper than that gas 
temperatures remain always CR-controlled. The small ~10% rise of T k (A v ) at the highest densities and Uqr values for A v >5 is likely an artifact of the 
high optical depths of C and CO at such large cloud depths (see Bayet et al. 2010). 



of high Ucr values on the initial conditions of star formation (Pa- 
padopoulos 2010). The thermal state of UV-shielded dense gas is 
indeed dramatically altered in the CRDRs of extreme starbursts, 
and this will affect its Jeans mass, gas fragmentation properties, 
and the charecteristic mass of young stars. This contradicts 
the results by Elemgreen et al. (2008) which, while correctly con- 
sidering the UV-shielded gas phase as the one where the SF initial 
conditions are set, it omitted CRs from its thermal balance. 

It must be pointed out that while several other studies investi- 
gated the effects of starburst environments on the IMF, all (besides 
the one by Elmegreen et al. 2008) failed to consider the ISM phase 
where the initial conditions of SF are actually set (the CRDRs). In- 
stead they used physical properties typical of PDRs (e.g. Klessen 
et al. 2007) which, while dominating the observed dust SEDs and 
molecular line Spectral Energy Distributions (SLEDs) of galaxies, 
they are hardly representative of the physical conditions expected 
in their CRDRs. Even for mildly SF IR-luminous galaxies the dust 
temperature typical of their global SEDs is ~(30-40) K, which is 



well above Td us t~(8-10) K expected in their CRDRs. Using inap- 
propriate PDR-type initial conditions of star formation in numerical 
models: a) would make them hard-pressed to account even for the 
near-invariant IMF found in the Galaxy (as the wide range of phys- 
ical conditions in its massive PDR gas phase would yield a wide 
range of IMFs), b) omits the significant, CR-induced, thermal de- 
coupling between gas and dust (Figures 2, 3). The latter fundamen- 
tally impacts the so-called inflection point of effective equation of 
state (EOS) used in all realistic (i.e non isothermal) numerical sim- 
ulations of molecular cloud fragmentation (e.g. Jappsen et al. 2005; 
Bonnell et al. 2006), and this is further discussed in section 3.2. 



3.1 The characteristic mass of young stars in CRDRs 

The thermal state of these dense pre-stellar gas cores is no longer a 
near-invariant but is instead strongly altered by the high CR energy 
density backgrounds of extreme starbursts, inevitably affecting the 
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Figure 3. The computed gas and dust temperature profiles versus optical extinction inside molecular clouds as in Figure 2, but now with the incident far-UV 
radiation field also scaling as the CR energy density. The now strongly varying dust temperatures are shown along with the gas temperatures as dashed lines 
with the same color. The gas temperatures at A v =5-10 remain almost invariant (compare with Figure 2) demonstrating the near-invariance of the thermal state 
of the gas (and thus of the star formation initial conditions) deep inside molecular clouds when it comes to its dependence on the radiation field. The abrupt 
Tfc drop at high densities (~(10 5 — 10 6 ) cm~ 3 ) and G values is due to a lack of atomic/ion coolants in the H-rich region: as H is converted into H2, the H2 
formation heating (which was the major heating agent) drops while at the same time CI cooling (CII is converted to CI due to H2/C mutual shielding) and H2 
cooling sharply increases. For high Ucr the gas remain much warmer than the dust, marking a CR-induced thermal decoupling of gas and dust (section 3.2). 



fragmentation of those gas regions towards protostars. Indeed, the 
Jeans mass for such CR-heated cores 

M (c) = / kBTjX 3/ %- 1/2 = 
V G/imH 2 J 





' n c (H 2 ) " 


-1/2 


V10K J 


_10 4 cm- 3 _ 





rises from M<y c) ~0.3M Q when n c (H 2 )=10 5 cm" 3 and T k =10K 
(typical for the Galaxy), to M^ c) ~(3-10) M for T k ~(50-1 10) K 
expected for same density H 2 gas immersed in the very large CR 
energy densities of extreme starbursts. The near-invariance of Mj 
in galaxies (Elmegreen et al. 2008) is thus upended in starbursts 
with high SFR densities. In Figure 4 the Mj values computed for 
molecular gas immersed in the CRDRs of compact extreme star- 
bursts are shown for a gas density range typical for star-forming 
molecular clouds. For UcR,/UcR,Gai^l0 3 , Mj increases by a fac- 
tor of ~10 across the entire range, which will invariably lead to 



a higher and thus a top-heavy IMF (Larson et al. 2005; 

Elmegreen et al. 2008). The wide range of densities over which this 
occurs makes this conclusion independent of the specific details 
of molecular cloud fragmentation as long as Mj [n(f, t), T k (f , t)] 
drives the fragmentation process outcome at each spatial and tem- 
poral point (f , t). The latter is the current consensus on the role of 
the Jeans mass in molecular cloud fragmentation towards a stellar 
IMF, though the views regarding the particular gas phase whose 
Mj sets the M^ vary (e.g. Klessen 2004; Bonnell et al. 2006). 

If M*>oc(Mj) at the onset of molecular cloud fragmentation 
(Klessen 2004), the CR-boost of in extreme CRDRs is ob- 
vious from the large vertical displacement of the Mj values in 
Figure 4. In a non-isothermal gravoturbulent molecular cloud frag- 
mentation scheme (e.g. Jappsen et al. 2005) it is M^'~Mj(n c ) at a 
characteristic gas density n c above which efficient gas-dust thermal 
coupling lowers the gas temperature in CRDRs towards a minimum 
value of T k (min)~Tdust (Larson 2005). For the Galaxy this occurs 
at ri c ~10 5 cm -3 , a density that marks also the transition to near- 
isothermal cloud regions where thermal motions dominate over su- 
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Figure 4. The Jeans mass estimated from Equation 3, and gas tempera- 
tures inside CRDRs, computed for the corresponding densities (using the 
full UCL-PDR code) and averaged over the inner A v =5-10 cloud regions 
(Figure 2). Low mass star formation will be clearly surpressed in extreme 
CRDRs for a Mj-driven gas fragmentation towards a stellar IMF (see text). 

personic turbulence (Larson 2005; Bergin & Tafalla 2007). From 
Figures 2, 3 it is obvious that the high CR energy densities in ex- 
treme CRDRs will keep the gas and dust thermally decoupled, with 
Tk(min)^>T dust, and thus a characteristic density can no longer be 
marked by the onset of thermal equilibrium between gas and dust. 
Such CR-induced effects have been long suspected for the dense 
molecular gas in the Galactic Center (Yusef-Zadeh et al. 2007). 

In the CR-innudated molecular clouds of extreme starbursts 
a characteristic density in the ISM can now remain only as the 
one where the supersonic turbulence has fully dissipated and over- 
taken by thermal motions in pre-stellar cores, a prerequisite for 
their efficient gravitational collapse (e.g. Krumholz & McKee 
2005). The average density of such transition regions can be esti- 
mated from two well-known scaling relations in molecular clouds: 
crv(r)=<7o(r/pc) h (linewidth-size) and (n)=n (r/pc) _1 (density- 
size) after setting <Tv(min)=(3kBTk//^mH 2 ) 1/ ' 2 as the minimum 
linewidth possible, and solving for the corresponding mean vol- 
ume density, 

<ntr)=n °U^r-) -MkrrTFr) {tok) ■ 

(4) 

For h=l/2 (expected for gas clouds virialiazed under ex- 
ternal pressure, Elemgreen 1989), o- =l.lkms" 1 , and for 
n o =fewxl0 3 cm" 3 , it is (n tr )~fewx 10 4 cm -3 . For the CR- 
boosted gas temperatures in extreme starbursts where typically 
T k (min)~100K (Figures 1, 2, 3): (n tr )~fewx 10 3 cm" 3 , and the 
(turbulent gas)— K thermal core) transition will occur at lower densi- 
ties as well as higher (CR-boosted) gas temperatures than in quies- 
cent CRDRs. This will then shift the expected M^*' both vertically 
(higher temperatures), as well as to the left (lower densities) of Fig- 
ure 4, yielding even higher values than those resulting from 
only a vertical shift of Mj because of CR-boosted gas temperatures. 

Thus irrespective of molecular cloud fragmentation details, a 
much larger characteristic mass of young stars is expected in the 
CR-innudated molecular clouds of extreme starbursts. However the 



exact shape and mass scale of the new stellar IMF can only be deter- 
mined with dedicated molecular cloud fragmentation simulations 
that make use of the new SF initial conditions in extreme CRDRs, 
along with an EOS appropriate for the very different thermal and 
ionization state of their molecular gas. 



3.2 CRDRs: SF initial conditions, the EOS inflection point 

The importance of the Effective Equation of State (EOS), 
parametrized as a polytrope V=Kp 1 , in the fragmentation proper- 
ties and emergent mass spectrum of turbulent self-gravitating gas 
has been well-documented in a few seminal papers (Li et al. 2003; 
Jappsen et al. 2005; Bonnell et al. 2006; Klessen et al. 2007) re- 
vealing a crucial role for 7 in defining the shape of the collapsed 
gas core mass distribution (and hence the stellar IMF). The EOS, 
its polytropic index 7 as a function of density, and their depen- 
dance on ISM properties such as metallicity, and physical condi- 
tions such as the elevated CMB at high redshifts have also been 
well-documented (Spaans & Silk 2000, 2005), with the so-called 
EOS inflection point emerging as a general characteristic. In an 
EOS that seems well-approximated by a piecewise polytrope with 
two distinct 7 values, the latter is simply the characteristic density 
n c for which y< 1 when n<n c , and 7>1 for n>n c . Its crucial role 
in determining the IMF mass scale, discussed briefly in 3. 1, is docu- 
mented throughout the literature (Larson 2005; Jappsen et al. 2005; 
Elmegreen 2008). In starburst environments a "sharpening" of that 
inflection point is thought to occur because of warm dust. There 
7 reaches up to ~1.4 and stays above unity for densities past the 
inflection point n c ~(10 4 -10 5 ) cm" 3 (Spaans & Silk 2005), which 
surpresses molecular cloud fragmentation while "tilting" it towards 
top heavy IMFs (e.g. Li et al. 2003; Klessen et al. 2007). 

Unfortunately PDR-like ISM has been used to set the SF 
initial conditions for starbursts in all of the aforementioned models, 
even though the dust in CRDRs (where these conditions are really 
set) remains cool, much lower than the gas, especially in extreme 
CRDRs. This insensitivity of the CRDR dust temperatures to the 
ambient far-UV radiation field can be demonstrated analytically 
by simply considering the location of CRDRs deep inside far-UV 
illuminated molecular clouds at A v > 5 (e.g. McKee 1989). Then 
Tdu8t,SB/Tdust,Gal~[eXp(-Av/1.086) x Go,sB/G ,Gai] 1/(4+a) 
will be the radiation induced boost of the CRDR dust tempera- 
ture in starbursts Td us t,SB with respect to that in the quiescent 
Galactic environments Td us t.Gai~8 K. For a dust emissivity 
law of a=2, a large enhancement of the far-UV radiation field 
Go,SB/G o ,Gai=10 4 , and Av=5 it is T dust ,sB/Tdust,Gai=2.15, 
yielding Td us t~17K, in good agreement with the detailed dust 
temperature profiles past A v ~5 shown in Figure 3. It is actually 
this robustness of dust temperatures inside CRDRs, along with a 
strong gas-dust coupling at high densities (when Ucr values are 
modest) that helps retain: the thermal states of the pre-stellar gas 
cores, the M^ values, and the IMF, as near invariants over a wide 
range of ambient PDR conditions (Elmegreen et al. 2008) 

Thus PDR conditions, while always dominating the observed 
dust SEDs and molecular SLEDs of even mildly SF galaxies (as 
PDRs are warm and contain the bulk of the ISM mass) they are not 
representative of the SF initial conditions in galaxies, which are set 
in CRDRs. This leaves the pre-stellar gas phase in them, amount- 
ing to a few% of the mass per Giant Molecular Cloud (~SF effi- 
ciency per typical GMC), yet to be explored in terms of its EOS 
and fragmentation properties. A detailed exploration of those in 
CRDRs will be the subject of a future paper, but Figures 2 and 
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Figure 5. The ratio of dynamical/cooling timescales in UV-shielded molec- 
ular gas of CRDRs (A v =5-I0) estimated using the UCL-PDR code and cos- 
mic ray energy densities ranging from those expected in quiescent Galactic 
to compact extreme starbursts. 



3 already show that CRs, by inducing a strong decoupling of the 
gas and dust temperatures with Tk3>Td us t, they "erase" any puta- 
tive EOS inflection point for gas densities up to n=10 6 cm~ 3 . This 
will maintain efficient molecular gas cooling, with cooling times 
always much shorter than dynamical times (Figure 5), and 7< 1 in 
CRDRs over the entire gas density range explored here. In such 
regions gas gravoturbulent fragmentation will then occur equally 
efficiently from n~5 x 10 3 cm -3 up to ~10 6 cm -3 yielding a pure 
power law IMF (Spaans & Silk 2000; Li et al. 2003). 

The latter result, along with our discussion in 3.1 makes clear 
that the stellar IMF will be very different in extreme CRDRs, 
though only detailed gravoturbulent simulations with the appropri- 
ate initial conditions, and performed with a new suitable EOS can 
decide its exact shape and mass scale. It must be also noted that the 
very high ionization fractions in extreme CRDRs can now "anchor" 
the magnetic field for much longer periods, and at much higher 
densities, as ambipolar diffusion in such regions is now much pro- 
longed (Papadopoulos 2010). Magnetic fields may thus remain im- 
portant during much later stages of gas fragmentation, even in a 
highly turbulent medium, and may thus have to be included in the 
simulations of this process in extreme CRDRs. 



4 CONSEQUENCES, OBSERVATIONAL PROSPECTS 

The current study supports the recently suggested role of CRs as a 
decisive SF feedback factor, capable of altering the SF initial con- 
ditions and the IMF in extreme CRDRs. In this section we discuss 
some of the consequences beyond the immediate issues of the ISM 
and the stellar IMF in such enviornments, and outline current and 
future observations that can shed light on these new SF initial con- 
ditions that will be typical in high SFR density events in the Uni- 
verse. 



4.1 A sequence of p s fr -dependent stellar IMFs 

The dependance of the fhermo-chemical and ionization state of pre- 
stellar dense gas in CRDRs on the average Ucr permeating the 
ISM of a galaxy naturally links the SF initial conditions and the re- 
sulting IMF to the average SFR density p s f r , provided Ucr oc p s f r . 
The details of CR transport/escape mechanisms in quiescent SF 
disks and starbursts will define the proportionality factor (see Suck- 
hov et al. 1993 for detailes). Thus, once simulations of molecular 
gas fragmentation in CRDRs produce libraries of the corresponding 
IMFs, these will be naturally parametrized by p s f r , and hence ac- 
tual observables (especially in the upcoming era of high-resolution 
imaging at submm wavelengths with ALMA). 



4.2 Implications for hierarchical galaxy formation models: a 
physically-motivated IMF bimodality 

In hierarchical galaxy formation models a bimodal IMF, top- 
heavy for mergers/starbursts and regular for disk/regular-SF sys- 
tems, has been postulated to explain the populations of the 
starbursting/merger-driven submillimeter galaxies (SMGs) and reg- 
ular SF disks like the Ly-break galaxies (LBGs) at high redshifts in 
ACDM cosmology (Baugh et al. 2005). The CR-controlled SF ini- 
tial conditions in galaxies, unhindered by dust extinction, can now 
set this on a firm physical basis. From Figures 2 and 4 is obvi- 
ous that for typical SF disks (UcR~(0.2-10)xUcR,Gai) no sub- 
stantial effects on the gas temperatures and Jeans masses in CR- 
DRs are expected while in the extreme CDRDs of compact star- 
bursts (Ucr~(10 3 -10 4 ) x Ucr, Gat) very warm gas and large Jeans 
masses will lead to an effective surpression of low-mass star forma- 
tion and hence a top-heavy IMF. For merger-driven star formation 
the latter may characterize entire star-forming episodes (i.e. until 
the in-situ molecular gas reservoir is nearly exausted) as dissipa- 
tion can continously "funnel" the large molecular gas masses of the 
progenitors towards very compact (-~(100-500) pc) regions where 
extreme, Eddington-limited SF, takes place with tremendeous IR 
brightness of eriR~(few)x 10 13 Lq kpc~ 2 (e.g. Thompson 2009). 
Such high IR brightness, and thus high SFR densities, have been 
recently deduced also in a distant SMG (Swinbank et al. 2010). 

Thus CR-controlled SF initial conditions naturally yield a bi- 
modal SF process in galaxies, with high p s f r values providing the 
branching point, as follows: 

(merger-driven SF)— Khigh p'sfr)— Khigh Ucr.)— Ktop heavy IMF), 
(regular disk SF)— >(low p s f r )— 5>(low Ucr)— ^(Galactic IMF). 
This proposed p s f r -IMF link can then parametrize such an SF 
bimodality in hierarchical galaxy formation, once IMF=F(p s f r ) 
"libraries" are obtained from e.g. gravoturbulent simulations of 
molecular cloud fragmentation for a range of p s f r values. 

A top-heavy IMF in high-p s f r systems will have very seri- 
ous implications on the actual SFRs deduced from the IR lumi- 
nosities of these typically very dust-enshrouded systems (where all 
the far-UV light from massive stars is reprocessed into IR). Their 
CR-induced top-heavy IMFs will now correspond to several times 
higher energy outputs per stellar mass, and the SFRs in such sys- 
tems will have to be revised downwards by identical factors. 



4.3 A p s f r -dependent IMF: integrated galactic IMFs 
dependent on star formation history? 

A p a f r -dependent IMF as the natural outcome of CR-controlled 
SF initial conditions in galaxies can also naturaly yield integrated 
galactic IMFs (IGIMFs) that depend on the star formation history 
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(SFH) of galaxies. This is because p slr can be a strong function 
of the evolution of a galaxy, especially during its early and pre- 
sumably gas-rich state even when no mergers are involved. Such 
a SFH-dependance of the IGIMFs has been considered before as 
the cause of the well-known mass-metallicity relation of galaxies 
(Koppen, Weidner, & Kroupa 2007). 

4.4 CR propagation in dense molecular gas: a dynamic view 
is needed 

Setting UcROcpstr in CRDRs assumes that CRs can freely pene- 
trate and warm the entire volume of a GMC down to its densest 
UV-shielded regions. Recent gamma-ray observations have directly 
measured CR energy densities in the nucleus of a nearby galaxy on 
par with those computed for extreme CRDRs and expected from its 
SFR density, while indicating only mildly calorimetric conditions 
for the CR/ISM p-p interactions (Acero et el. 2009). The latter sup- 
ports an unhintered CR penetration through the molecular clouds. 
Nevertheless Alfven waves generated outside such clouds by a net 
flux of CRs in their environments could be effective in keeping 
low energy CRs below a few hundred MeV (those most effective 
in heating the gas) outside their high density regions (Skilling & 
Strong 1976), which in turn would affect their thermal balance, 
their fragmentation properties, and the emergent stellar IMF. Thus 
the coupled problems of CR transport, magnetic field dynamics, 
and gravoturbulent evolution of non-isothermal dense molecular 
gas with high degrees of ionization like those expected in compact 
extreme starbursts such as Arp 220 (e.g. Greve et al. 2009) must be 
investigated with dedicated MHD simulations. 

4.5 Observational prospects 

The detections of 7-ray emission from two nearby galaxies, namely 
M 82 and the starburst nucleus of NGC 253 (Acciari et al. 2009, 
Acero et al. 2009), are two very important recent developments in 
7-ray astronomy, given that 7-ray emission results from CR proton 
interactions with the molecular ISM. This allows such observations 
to directly measure the average CR energy densities in SF environ- 
ments, and thus sensitive such measurements in ULIRGs such as 
the nearby ULIRG Arp 220 will be very important in revealing ex- 
treme CRDRs in such systems. 

On the other hand, direct measurements of the expected top- 
heavy IMFs in the heavily dust-enshrouded compact starbursts 
where extreme CRDRs would reside are impossible. However any 
cases where the "end" products, i.e. the stellar populations, of high- 
psh SF events can be observed and their IMF determined such as 
the Galactic Center (GC) and local Ultra Compact Dwarfs (UCDs) 
allow such direct observational tests. There is indeed strong evi- 
dence of CR heating of the dense molecular gas in the GC, setting 
Tki n >Tdust, and expected strong effects on the local Jeans mass 
and ambipolar diffusion timescales (Yusef-Zadeh et al. 2007). Un- 
fortunately, initial assertions about a top-heavy IMF in the GC (e.g. 
Nayakshin & Sunyaev 2005) and the massive Arches cluster near 
it (Figer et al. 1999; Stolte et al. 2002, 2005) have not been cor- 
roborated by later studies (Kim et al. 2006; Espinoza et al. 2009; 
Lockmann et al. 2010). On the other hand some evidence exists for 
a top-heavy IMF in the nuclear starburst of M 82 (e.g. Rieke et al. 
1980; 1993), and a significantly enhanced Ucr is certainly there 
(~500xUcr,Gb1, directly measured via 7-ray observations; Ac- 
ciari et al. 2009) to achieve this via a CR-boosted IMF mass 
scale. Indeed for such Ucr values, setting M^'~M.j(10 5 cm~ 3 ) 



yields M);* , (M82)>2Mo (Fig.4). However the Galactic Center 
does offer an example of regions with locally boosted Ucr val- 
ues that do not necessarily correspond to top-heavy IMFs in their 
vicinity, possibly because such boosts are too short and transient to 
systematically bias the local emergent IMF. 

Stronger evidence for top-heavy IMFs exists in UCDs, 
small stellar systems with r e ~(10-30)pc (half-light radius) 
where ~(10 6 -10 7 )M© of gas have been converted into stars 
at (SFR)~(lO-lOO)M yr _1 (see Dabringhausen et al. 2009 
and references therein). Unlike the GC where IMF studies are 
hindered by high extinction and tidal disruption of clusters, 
and where any high-p slr SF episodes may have been short- 
lived (and thus could not significantly bias its average IMF to- 
wards a top-heavy one), the IMFs in UCDs are much more 
amenable to systematic studies, while their SF episodes have had 
(p sfr (UCDs))>(p Bfr (ULIRGs))! They are thus ideal targets for 
studying the IMF resulting from completed high-p slr SF events. 

The next best thing besides studies of the IMF itself is stud- 
ies of the ISM in CRDRs from which the IMF emerges, especially 
those in compact extreme starbursts where p slr is very high (e.g. 
ULIRGs) versus systems with more distributed SF and much lower 
p s f r values. The unique observational signatures of the very warm 
dense gas with high ionization fractions expected in extreme CR- 
DRs have already been discussed in some detail by Papadopoulos 
(2010), and the interested reader is refered to this work. In sum- 
mary, in the era of Herschel and the upcoming era of ALMA the 
sole main obstacle is the degeneracy of extreme CRDR diagnos- 
tic with that of X-ray Dominated regions (XDRs) induced in the 
ISM by an X-ray luminous AGN. Indeed, unless an X-ray luminous 
AGN can be excluded by other independent means (e.g. deep X-ray 
observations), it can be difficult to tell extreme CRDRs and XRDs 
apart in a simple fashion as even exotic high excitation species such 
as OH + and H20 + , recently detected by Herschel in the archety- 
pal ULIRG/QSO Mrk231 and attributed to XDRs (van der Werf 
et al. 2010), could be partly due to CRDRs. In recent work Mei- 
jerink et al. 2010 provide a detailed discussion of diagnostics lines 
in extreme CRDRs versus XDRs and the key observations that can 
separate them. 

Finally it is worth noting that X-rays and XDRs will have sim- 
ilar effects as CRs on the Jeans mass and the IMF of UV-shielded 
dense gas regions as both can volumetrically warm such regions to 
high temperatures. This was pointed out recently by Bradford et al. 
(2009) for the ISM of a distant QSO, and constitutes an omitted 
aspect of AGN feedback on star formation. This has been studied 
recently and shown to lead to a top-heavy IMF (Hocuk & Spaans 
2010). XDRs besides elevating Tk(min) to much higher values than 
in the Galaxy they also yield Tki n >-Tdust (and thus they can also 
"erase" the EOS inflection point), though unlike CRDRs their influ- 
ence on the ISM is necessarily limited by the ~1/Ragn geometric 
factor (Ragn the distance of the irradiated gas from the AGN). 

5 CONCLUSIONS 

We have conducted new, detailed calculations on the thermal bal- 
ance of UV-shielded dense gas in Cosmic-Ray-Dominated-Regions 
(CRDRs), the ISM phase where the initial conditions of star forma- 
tion are set in galaxies, in order to examine in detail the effects of 
extreme CRDRs on the mass scale of young stars and the IMF re- 
cently suggested by Papadopoulos (2010). Our results can be sum- 
marized as follows: 

1 . We confirm the high temperatures of dense gas cores in the 
extreme CRDRs expected in ULIRGs and all systems with large 
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star formation rate densities (p s f r ), albeit at lower levels than the 
original study. 

2. These CR-induced higher gas temperatures will lead to ~10 
times larger Jeans masses across the entire density range typical in 
such regions, and thus higher charecteristic mass of young stars in 
the ISM of galaxies with high p s f r . 

3. The CRs, by decoupling the thermal state of the gas from 
that of the dust in CRDRs keep Tk3>Td us t for all densities typical 
of molecular clouds. This effectively "erases" the so called inflec- 
tion point of the Effective Equation of State, especially in extreme 
CRDRs, inevitably altering the fragmentation properties of the gas. 
This effect, along with the expected much higher gas ionization 
fractions, necessitates a search for a new EOS, appropriate for CR- 
DRs. 

4. These new CR-induced ISM conditions provide a new 
paradigm for all high-density star formation in the Universe, with 
CRs as the key SF feedback mechanism, operating unhindered by 
extinction even in the most dust-enshrouded starbursts. 

5. The resulting p s f r -dependent IMF, yields: a) a natural bi- 
modal behaviour for merger-driven/starburst versus disk/regular-SF 
galaxies, and b) integrated galactic IMFs (IGIMFs) that depend on 
their star formation history, both with important consequences in 
galaxy formation and evolution. 

6. IMF studies of stellar systems formed under extremely high 
(p s fr) values such as Ultra Compact Dwarfs, and molecular line ob- 
servations of the extreme CRDRs in nearby ULIRGs can provide 
the full picture of CRDR properties and their effects on the IMF. 
Regarding the latter, the now spaceborn Herschel Space Observa- 
tory, and the upcoming ALMA hold particular promise. 
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APPENDIX A: THE COOLING AND HEATING FUNCTIONS FOR UV-SHIELDED CORES 

The CR heating rate used in this work is 

Tor = 1.95 x lCr 28 n H ( — < ^— 7 r ) crgscm-^- 1 (Al) 

\ 1.3 x 10~ 17 s _1 J 

(Wolfire et al. 1995) where CcrocUcr being the Cosmic Ray ionization rate (in s -1 ), with a adopted Galactic value of (xir=5x10 -17 
s _1 (corresponding to UcR=UcR,Gai), and nH=2n(H2) for fully molecular gas cores. For an optically thin OI line the cooling is Aoi63 = 
XoiihCiuE u i, which becomes 

Aoiea = 3.14 x lCT 14 xon H C ul {^j exp (-227.72/T k ) ergs (A2) 

where g u =3 and gi=5, and xo=[0/H] being the abundance of oxygen not locked onto CO (xo~4.89x 1CP 4 ). The collisional de-excitation 
coefficient is given by 

C ul = nHl0°' 321ogTk - 10 ' 52 = 3.02 x 10- 11 n H 'lt- S2 cm- 3 B- 1 , (A3) 
(Liseau et al. 1999). Thus finally Equation A2 becomes, 

Aoiea = 2.78 x 10~ 28 nHT°' 32 exp (-227.72/T k ) ergs cm^s" 1 . (A4) 
At lower densities (nu <10 5 cm -3 ) and strong CR fluxes, a small fraction of carbon remains in the form of CII, acting as a coolant with 

Aon = 1.975 x l(T 23 n!xcii exp (-92.2/T k ) ergs cm" 3 s _1 , (A5) 

computed in a similar fashion as the OI(63 fim) line cooling, and for a fully transparent medium. Gas-grain accommodation can heat or cool 
the gas depending on the difference between their temperatures (Burke & Hollenbach 1983). It can be expressed as 



T acc = 4.0 x 10 12 an H n d cr g ^ (T k - T dust ) (A6) 

where n g is the number density of dust grains and er g is the grain cross-section (cm -2 ). If the gas-to-dust ratio is 100, the dust mass density 
is pd=3.5 g cm~ 3 , and a dust radius of a — 0. 1 fim then 

n d = 2.2 x 10~ 2 x ^ H n H /(4/37rp da 3 ) = 7.88 x 10~ 12 n H . (A7) 
For the accommodation factor a we follow the treatment by Groenewegen (1994) where 



a = 0.35 exp (^- ] fl^±I^, (A8) 

which we set to the maximum value a=0.35 (i.e. maximum gas cooling from gas-dust interaction). Thus the mainly cooling term (for the 
ISM conditions explored here) due to the gas-dust interaction in Equation A6 becomes 

A gd = 3.47 x 10 _33 n^A/Tk (T k - T dust ) ergscnT 3 s" 1 . (A9) 
Finally based on our UCL-PDR code outputs we parametrize the cooling due to the CO rotational transitions as 



where we set Xco/X[C] =(0.97,0. 98,0.99, 1.0) for densities of 5xl0 3 , 10 4 , 10 5 , 10 6 cm -3 respectively (with the value of =1 corresponding 
to all carbon locked onto CO). 

Solutions of Equation 2 in the main text for densities spanning the range present in star-forming H2 clouds in galaxies, and Ucr=(0.2, 
1.0, 10, 10 2 , 10 3 , 10 4 )xUcR,Gai are then obtained numerically and are shown in Figure 1. 



